Application of Bryan's algorithm to the mobility spectrum analysis of semiconductor 

devices 
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A powerful method for mobility spectrum analysis is presented, based on Bryan's maximum 
entropy algorithm. The Bayesian analysis central to Bryan's algorithm ensures that we avoid over- 
fitting of data, resulting in a physically reasonable solution. The algorithm is fast, and allows the 
analysis of large quantities of data, removing the bias of data selection inherent in all previous 
techniques. Existing mobility spectrum analysis systems are reviewed, and the performance of the 
Bryan's Algorithm Mobility Spectrum (BAMS) approach is demonstrated using synthetic data sets. 
Analysis of experimental data is briefly discussed. We find that BAMS performs well compared to 
existing mobility spectrum methods. 

PACS numbers: 72.20.My, 02.70.-c, 95.75.Pq 



I. INTRODUCTION 

The characterization of semiconductor devices is es- 
sential for the development of new materials which 
are interesting from both a physical and technological 
point of view. Hall-effect measurements are very fre- 
quently used for this characterization, especially when 
two-dimensional systems are investigated. 1] 

The simplest ideal material has a single charge car- 
rier with isotropic energy bands and energy-independent 
scattering mechanisms. The resulting longitudinal resis- 
tivity p and Hall coefficient Rh are constant with respect 
to the application of magnetic field. Such a picture is 
effective in a basic understanding of certain solid state 
phenomena. 

In modern semiconductor devices, heterostructures are 
increasingly employed where different semiconductor ma- 
terials are grown epitaxially in layers on a substrate. In 
this manner, structures can be designed in which charge 
carriers are confined in a quantum well physically remote 
from the dopant atoms they originate from; ionized im- 
purity atoms can strongly limit the mobility of charge 
carriers so this method of modulation dovingc&n lead to 
greatly enhanced transport properties. 0,0)13 

However, the physics of magnetotransport in such het- 
erostructures is considerably more complicated if there 
are two or more distinct carrier gases present in the mate- 
rial (for example, the intended carrier gas in the quantum 
well plus the doped region) or the carrier gases feature a 
spread of mobilities due to energy-dependent scattering 
mechanisms or multiple subband occupancy. This leads 
to resistivities and Hall coefficients which are dependent 
on the applied magnetic field, and to extract the prop- 
erties of such systems the simple single-carrier model is 
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not sufficient. To this end, various mobility spectrum 
techniques have been employed. 0, 0, 0, El 

In a seminal paper, based on the work of McClure,0 
Beck and AndersonQ showed that the conductivity ten- 
sor depends on a generalized conductivity function, s(p,), 
via the integral transforms, 



and 
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and n s (p) represents the number (per unit area) of car- 
riers with a mobility p averaged over the sample. 

The Hall coefficient and resistivity are related to the 
elements of the conductivity tensor in the standard way, 



p{B) = pxx 



o~xx 



'xx 



Rh{B) = - 



PXY 

B 



' XY 



0~XY 



(4) 



(5) 



'XY 



The function s(p) appearing in Eqs. Q-© is often 
known as the mobility spectrum. It gives the contribu- 
tion to the conductivity due to the density of carriers 
with mobility p. In this investigation, it is assumed that 
s(p) does not change with magnetic field. At high fields 
and low temperatures this is not true because of the for- 
mation of Landau levels: the density of states becomes a 
function of magnetic field if Huj c > k^T where oj c is the 
cyclotron frequency eB/m*. 9\ Also, at very low fields 
and low temperatures (such that pB < h/Epr where 
Ep is the Fermi energy and t is the momentum relax- 
ation time) quantum effects such as weak localization 
contribute to the magnetoresistance. [Toj 
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According to Eqns. and J5J magnetoresistance is 
always positive; negative magnetoresistance can only be 
analyzed by assuming a magnetic-field dependent mobil- 
ity or momentum relaxation time. [Til Il2| 

Finding the magnetoresistance given a form for the mo- 
bility spectrum is straightforward. However, performing 
the so-called inverse transform problem (that is, finding 
the mobility spectrum from magnetoresistance data) is 
non-trivial with a noisy, limited data set and is the sub- 
ject of much investigation. @, EH E3, EM EH E3, E^l 

We replace the traditional MAXENT algorithm used 
by Kiatgamolchai et al. 0, 0] with Bryan's 
algorithm. 20] This approach has two major advantages 
over the previous technique: (1) A Bayesian analysis is 
performed on the resulting spectrum, and (2) the search 
process is optimized by null space decomposition of the 
kernel. The Bayesian aspect of Bryan's algorithm ensures 
that we choose the most probable balance between the 
fit quality, given by least squared minimization, and the 
information content of the solution, given by the maxi- 
mization of entropy. The null space decomposition pro- 
vides an exponential basis for the solution. The expo- 
nential basis is particularly well suited to this type of 
problem since it significantly reduces the number of di- 
mensions in the search space. It therefore allows reliable 
well-conditioned analysis of greater quantities of data, 
avoiding the need for potentially biased data selection. In 
this way, the results should be relatively free from com- 
mon artifacts such as mirror peaks, peaks of unnatural 
shape, or structure in the spectrum at very low mobility 
values, provided that the artifacts are not introduced by 
systematic measurement errors, pi) 

We present this paper as follows. In Sec. [HJ early ap- 
proaches to this problem are reviewed, before consider- 
ing the method of maximum entropy in more detail. In 
Sec. IIIII the maximum entropy method used in this paper 
is introduced, and technical details of Bryan's algorithm 
are described. In Sec. IIVI we perform a case study on 
model data using our Bryan's Algorithm Mobility Spec- 
trum (BAMS) method. The usefulness of the mobility 
spectrum technique for the analysis of real data is dis- 
cussed in light of these results, and a brief summary is 
given in Sec. [3 



II. BACKGROUND 

The approach of Beck and Anderson, Q which first 
introduced the concept of the mobility spectrum and 
has been developed into commercial software, involves 
a multi-stage solution method (the numerical details of 
which will not be explored). The number of data points 
used to create a solution is equal to the number of dis- 
tinct carriers gases within the sample, plus one. The 
method either searches for an imposed number of car- 
rier gases or fits for a single carrier system, then a two- 
carrier system, and so on. During each search for n gases, 
the method tries fitting to every available combination of 



n + 1 data points from the complete data set, only keep- 
ing the sets which lead to "physical" solutions. From 
the retained sets, the algorithm generates an "envelope 
function" within which the mobility spectrum should be 
contained. 

Once a set of carrier gases has been found, the validity 
of the solution is confirmed by performing a least-squares 
fit on the original data with respect to the carrier gas pa- 
rameters. If the parameters emerge virtually unchanged 
from the second stage, then the results of the original 
analysis are considered to be reasonable. 

If only a limited number of data points (~ 10) are avail- 
able and the system contains only two or three distinct 
(i.e. with mobilities which differ significantly from each 
other) carrier gases, then this method is quite efficient. 
However, with automated systems the data collected can 
comprise thousands of values and the system simply can- 
not cope; choosing points or creating averaged data by 
hand imposes artificial constraints on the final solution. 
Also, this method does not respond well to errors in the 
experimental data.fl8| 

Reference ^| introduces and describes the Reduced 
Conductivity Tensor method for extracting the carrier 
concentration and mobility of each component of a mul- 
tilayer semiconductor system. It is assumed that the car- 
rier gases are essentially degenerate and that the effec- 
tive mass is isotropic. These assumptions immediately 
remove the possibility of gaining insight into energy- 
dependent scattering mechanisms in non-degenerate sys- 
tems, and also the possibility of analyzing for exam- 
ple p-type conduction in pure germanium in which the 
anisotropy of the heavy hole band should be evident as 
"harmonics" in the mobility spectrum. 0, I22I l23j Nev- 
ertheless, the method works well for up to three carrier 
gases. 

In Ref. 15 it is assumed that the mobility for a gas 
subjected to mixed scattering can be approximated by 
the phenomenological expression: 

E 

H = fi x a where x = —— (6) 
k B T 

and that the mobility spectrum of non-degenerate carri- 
ers in a spherical band can be approximated by: 

s(M) = Soe-'Mxfr) (7) 

where So, /xo and a are coefficients depending on the 
density of states and on the parameters of the scatter- 
ing mechanisms. Then, the integrations in Eqs. Q and 
P)l are converted to summations of discrete spectra for 
electrons and holes and an iterative transformation pro- 
cedure may be performed. The results are decomposed 
into high and low-mobility carrier contributions. The dis- 
tinction between low-mobility carriers and high-mobility 
carriers is to some extent arbitrary, however, it corre- 
sponds to the visible separation between two different 
regions on the mobility spectrum. Results for the high- 
mobility carrier gas seem much more satisfactory than 
results for the low-mobility carrier gas in each case. 
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References and ^3 describe the quantitative mobil- 
ity spectrum technique (QMSA) and the improved quan- 
titative mobility spectrum technique (i-QMSA) respec- 
tively. These are also iterative techniques but with 
no initial assumption about the solution (although the 
Beck- Anderson Q approach is used to create a trial so- 
lution, and conductivity data is extrapolated to higher 
magnetic field values to extend the available mobility 
range). The i-QMSA method introduces a few extra 
tricks for improving the fits whilst smoothing the spec- 
tra and making them more physically reasonable. This 
method has also been successfully applied to conduc- 
tion in anisotropic bands by in corp orating an explicit 
anisotropy coefficient. [24j In Ref.llfll it is shown that this 
method does not deal with error very well. 



Pi = ^ (9) 
Co 

where pj is chosen so that the spectrum has a normaliza- 
tion of unity and <7o is the conductivity at zero magnetic 
field. The form of the entropy imposes the condition that 
the mobility spectrum is non-negative, which is physi- 
cally realistic. If any prior information about the form of 
the spectrum is available then this may be incorporated 
as a so-called default model Iml.poj 

Fitting is now a matter of minimizing the function 




Q = x 2 - aS (10) 



III. MAXIMUM ENTROPY 

A serious issue (regarding inverse transformation prob- 
lems in general) is that while a particular spectrum s(/i) 
may easily be transformed to produce a magnetoconduc- 
tivity cr(-B), the inverse problem is ill-conditioned and 
the solution obtained for s(/z) by inverting the kernel is 
extremely sensitive to small changes or errors in the o~(B) 
data. In particular, there is not necessarily a unique s(fi) 
within the uncertainty bounds of the original data. For 
example, since the integral transform has a blurring ef- 
fect, then it is possible for the data to have large local 
fluctuations while producing a smooth cr(B). 

The technique of Maximum Entropy has been em- 
ployed to solve this and other problems involving inverse 
integral transformation where the result is a positive, ad- 
ditive function. [H Him mill Generally, if two solu- 
tions (found by any means) are of equal merit (in terms 
of, for example, their least-squares fits to the original re- 
sults) then the solution with the larger "entropy" is to 
be favoured, since it is maximally non-committal with 
regard to missing (unmeasured) information in the origi- 
nal data. In other words, the solution favoured by 
the maximum entropy method extracts the most infor- 
mation out of the original data with the most reasonable 
assumptions regarding information which is unavailable. 

By discretising the mobility spectrum {s}, the entropy 
S can be defined as:|2i| 

N 
4=1 

I 



where \ 2 — \ J2i f ^i^-^^&l ^ ■ a c (B) is magnetocon- 
ductivity calculated from the fitting mobility spectrum, S 
is the error in the data such that a chi-squared which re- 
spects the data has a value of approximately M (the num- 
ber of "observations" or data points) and a is a hyper- 
parameter which controls the relative importance of the 
least-squares and entropic constraints. [20] 

There are various flavors of the maximum entropy 
method, which differ in the way a is determined. In 
the case of historic maximum entropy, different values of 
a are tried, until the minimization of Q [Eq. (|10|l ] gives 
X 2 = M. In the case of classic maximum entropy, the 
most probable value of a is found, given the data and 
the default model. [2^, |3(j Alternatively, the method de- 
scribed in Ref. treats the maximization of the entropy 
as a natural starting-point for derivations of other results 
of statistical mechanics. 

The approach of Bryan's algorithm is to calculate spec- 
tra for a range of a values and evaluate the probability 
that the solution is correct given the data and the default 
model. A weighted average over a is then taken. [20] How- 
ever, this method may not work well if a is very small. 
In the method developed in Ref. |3jJ a solution produced 
by Bryan's algorithm maximum entropy is used as the 
default model for the next solution. This procedure is it- 
erated until the most probable solution corresponds to an 
a value well within the range of applicability of Bryan's 
algorithm. 

Since there are two integral transforms involved in cal- 
culating cr(B), then we work with|l9j| 



1 1 f+^ 2 

a{B) = axx(B) + n XY {B) = -^o + 7: / s(/i)[cos(20) + sin(20)]d0 (11) 

1 1 J-tt/2 



where tan(0) = fiB. This produces a fit to uxx{B) + o~xy(B) rather than optimally to each component, which 
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could be seen as a drawback. Efforts towards the develop- 
ment of a method which inverts each of the o~xx(B) and 
<jxy{B) tranforms in a semi- independent manner have 
not been successful so far. However, since <Jxx{B) is an 
even function and o~xy{B) is odd, then if both negative 
and positive magnetic field values are used symmetrically, 
a fit to <Jxx{B) + <jxy{B) unambiguously fits <txx{B) 
and <jxy{B). This improves results, especially for low 
mobility carrier gases, as will be seen in Sec. IIVI 

Since we have assumed discrete forms for s(/i) and 
o~(B), then Eq. Ijllfl can be rewritten o~i = KijSj (with 
M data points and N mobility spectrum points). The 
mobility spectrum may be computed as: 



Minimising Q leads to the following relation, 



(12) 



Since the magnetoconductivity o~{Bi) is relatively in- 
sensitive to the details of the mobility spectrum s(fij), 
the kernel contains a large quantity of repeated infor- 
mation; many of the linear equations described by the 
kernel are (almost) identical. |33] In addition, elements of 
the kernel matrix which are dominated by signal noise or 
computational rounding errors are liable to make a huge 
contribition to Eq. (jl2|l . 

The dominant effects of data error in the problem can 
be investigated via singular value decomposition (SVD). 
The M x N kernel may generally be rewritten as K — 
UWV T where U is an M x N column-orthogonal matrix, 
W an N x N diagonal matrix with nonnegative values 
(the "singular values") and V T the transpose of an or- 
thogonal N x N matrix V. As a result of the orthogo- 
nalities, the inverse of the kernel can simply be written 
as K- 1 = VW~ 1 U T .^ 

If one of the elements of W is zero then the correspond- 
ing element of W~ l will be zero and the matrix is singu- 
lar. However, small (non-zero) elements resulting from 
inaccuracies invert to large values, whereas they should 
become zeros. It is therefore justified to set values of W 
with magnitude smaller than the data error to zero.|32j 
SVD reduces the dimensionality of the space that must 
be searched for a solution giving a massive increase in 
speed, and discards the effects of noise which would cor- 
rupt the solution. It is therefore a very important feature 
of the following analysis. 

A simple, powerful and successful maximum-entropy 
method for finding a mobility spectrum is described in 
Refs. [Hand [H based on the method of Javnes.[28j The 
described procedure does not, however, use SVD to re- 
move overspecification and noise and so the calculation 
time increases with the square of the product of the num- 
ber of mobility and magnetic field points, (A^M) 2 ;E3 nor 
does it make explicit a hyper-parameter [a in Eq. pOfl] 
for balancing fit quality against the maximization of the 
entropy, or deal with error in the data in an obvious way. 
Also, there does not appear to be any way in which the 
likelihood of the solution is checked and this leads to 
overfitting and the creation of artifacts in the mobility 
spectrum when applied to real data.[4| 
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dx 2 



aln(— ) - / t 
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(13) 



By introducing the SVD to this relation, it can be 
shown that the Sj can be cast in the form of an expo- 
nential basis, 



Si = mi exp(£ Uffiij) 

3=1 



(14) 



where the superscript (s) represents matrices in the re- 
duced space. 

The solution may then be found in terms of the s- 
dimensional Ut space rather than the 7V output -dimensional 
Si space. The Ut are found iteratively by using a Newton 
method to obtain the increment in parameter space Suj, 

^(a + e)8ij + ^2 M ik D kjJ 5uj = -aUj - gj (15) 



where, 



1 .»V^T/W r 9x2 



2~ ! ^ 13 da 3 {B) 



(16) 



and two new matrices, 



M = iw< s ) [V^f^V^W^ (17) 



D = [U (s) ] T diag[s]U 



(18) 



are defined to speed up the algorithm: M can be calcu- 
lated once at the beginning of the algorithm, leaving the 
smaller computation of D to the iterative process. [33j 

Many solutions are found for different values of a and, 
by considering the curvature of the parameter space in 
the vicinity of the solution, the relative probability of 
the spectrum for a given a can be found. Performing 
the weighted average should then give the most probable 
mobility spectrum. 



P[a[m, CI \ J ' 



a + Aj 



1/2 



exp[- X 2 /2 + aS]P[a] 



Aj are eigenvalues of the matrix, Aj , which describes the 
curvature of parameter space in the vicinity of the solu- 
tion, 



A = diagts^^MUdiagls 1 



/2l 



(20) 



and P [a] (which must be rescalable) is normally given by 
the Jeffries prior .pp] 



P[a] = 1/a 



(21) 
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In the BAMS method described below, the calculation 
time scales as N^M - 6 . The SVD scales as N 2 M, and 
remains a few percent of the calculation time. For input 
data at 500 magnetic field values and an output spec- 
trum of 300 points (finding solutions at 1000 values of a) 
the first BAMS result is produced in roughly 3 minutes 
on a 1.7GHz Intel Pentium 4 system running Linux, the 
first few seconds of which are the SVD. Such a number of 
input points is realistic given automated data-gathering 
systems; the time-limiting factor in such experiments is 
often the rate at which the magnetic field can be swept. 
All previous mobility spectrum analysis methods require 
that ~ 90 % of a typical data set were discarded by the 
user in order that the algorithm would complete in a rea- 
sonable time or indeed work at all. The SVD essentially 
performs this same task (of reducing the space that must 
be searched for a solution) in a mathematically rigorous 
manner, removing user bias and also saving time. 



1.0 - 
0.8 - 



£ 0.6 —| 

E 



3 0.4 —| 

03 



0.2 - 
0.0 



n 1 1 r 



T 

■ Actual 



t 1 1 r 



0.0 




u [ m VV ] 



1.0 



IV. RESULTS ON SYNTHETIC DATA 

In order to characterize the BAMS solution method, 
synthetic data were generated using Eqs. with re- 

alistic parameters and subject to varying degrees of noise. 
Carrier gases were modelled as Gaussian peaks in s(/i) 
centered at the specified values. Since the conductivity 
of the carrier gas is given by the area under the peak [see 
Eq. (|22|) ]. the intended sheet carrier density and peak 
width were also specified for each carrier gas. The height 
of a peak in s(n) generally has no physical significance. 

There were 1000 (B, (Txx,vxy) points from — 10 T to 
10 T in each case, apart from Fig. [3 where a dataset with 
500 points from T to 10 T is used for comparison. 

The results shown in Fig. ^ (f° r synthetic data 
created from a spectrum of two Gaussian peaks 
with fii=0.2 m 2 V _x s -1 and n x =l x 10 11 cm" 2 , 
/i2=0.6 m V s and n.2=l x 10 11 cm~ 2 ; in this case 
the width of each peak is proportional to the peak mo- 
bility) demonstrate how error in the data leads to very 
little broadening and inaccuracy of peaks in the output 
spectrum. The iterative procedure where one solution 
is used as the default model for the next [see Eq. JSJ] 
was not necessary in this case. By integrating over each 
peak (nominating a suitable [i v&a k value) the number of 
carriers can be found within each gas: 



peak 



s(/i)d/i = n s ejj, peak 



(22) 



Additionally, since the Hall sheet density is defined as 
nHaii = {eRn{B — * 0)) -1 then from Eqs. © and ©: 



«Hall 
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'.J jus(»d/x 



MHaii = <toRh{B -> 0) = — / /j,s(n)dn 

CO J 



(23) 



(24) 



FIG. 1: Mobility spectra calculated on 1000-point synthetic 
data featuring two carrier gases (the solid curve designated 
"Actual" is the specified spectrum) with the specified level 
of random data error. There were 1000 output points in 
the mobility spectrum between /i = -2.0 m 2 V -1 s -1 and 
+2.0 m 2 V — l s — 1 but only the region of interest is shown. Even 
3% error in the data does not spoil the result. Peak positions 
are summarized in Table [I] 



The true (drift) carrier concentration is simply found 
by summing the carrier concentrations of each peak: 



^Drift 



E 



so the drift mobility is: 



MDrift 



00 



(25) 



(26) 



A summary of the results from Fig.^are shown in Ta- 
ble [U a simple algorithm was used to find local maxima 
in the spectra. (Rather than fit a Gaussian to the spec- 
trum to determine the position of the peak, we simply 
find the highest point: this method is therefore limited 
in its precision by the number of mobility points chosen 
for the spectrum.) In addition to the peaks described, 
there were generally a few much smaller peaks contribut- 
ing in total to less than 0.5% of the conductivity. (For 
clarity, not all of the results summarized in Table [I] are 
shown in Fig. [I]) 

Generally, when the peak mobility is overestimated the 
sheet density is underestimated (and vice-versa) lead- 
ing to a correct partial conductivity. The Hall mobility 
compares well with the "ideal" value of 5000cm 2 V _1 s _1 
at 1.6xl0 n cm~ 2 (given by the two-carrier Hall-effect 
equations). Even with 3% error in the data the results 
are acceptable. 

With this software (as with that described in Refs. Il8l 
and 119) it is possible to find peaks in the mobility spec- 
trum at mobilities smaller than 1/B max , and this is 



6 



TABLE I: Summary of results of BAMS analysis on the 2-peak synthetic data shown in Fig.0 It can be seen that even up to 
a 3% data error the carrier gases are located accurately. Since the algorithm for finding the peak position is very simple, the 
precision is dependent on the number of mobility points. 
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FIG. 2: Mobility spectra calculated on synthetic data fea- 
turing two carrier gases of low mobility, with 0.01% random 
data error. There were 200 output points in the mobility 
spectrum between fi = -0.5 m 2 V -1 s -1 and +0.5 m 2 V _1 s _1 . 
It can be seen that the result changes depending on whether 
only one field direction (500 points) or both (1000 points) are 
used even though this formally does not change the amount 
of information. Peak positions are summarized in Table ITU 



shown in Fig. [21 and Table |nj The synthetic data fea- 
ture two carrier gases (^i=200 cm 2 V — 1 s _1 and n%=l x 
10 13 cm- 2 , ^2=1500 cm 2 V" 1 s- 1 and n 2 =l x 10 12 cm- 2 ) 
with 0.01% random data error. 

If only positive field values are used, the peaks are 
found at mobility (density) values of 377 cm 2 V _1 s _1 
(5.6 x 10 12 cm- 2 ) and 1533 cm 2 V -1 s _1 (9.0 x 10 11 cm" 2 ). 
The low mobility peak is not in exactly the correct posi- 
tion, but the high mobility peak (which tends to be more 
important for real systems) is located accurately. 

However, if both positive and negative field values 
are used the peaks are found to be much closer to 



the specified values, at mobility (density) values of 
226 cm 2 V- 1 s- 1 (8.8 x 10 12 cm" 2 ) and 1482 cm^^s- 1 
(1.0 x 10 12 cm -2 ). This is because the use of positive and 
negative field values fits better to axx and <jxy whilst 
the the use of one field direction only fits to oxx + oxy ■ 

In Fig. |31 which features three carrier gases, the evo- 
lution of the spectrum with increasing iteration over de- 
fault models is shown. The positions of the peaks are 
summarized in Table ITTT1 With increasing iteration the 
peaks become too narrow, but the accuracy of the de- 
termination of the peak position increases slightly. This 
demonstrates that many iterations over successive default 
models is not necessary; it also shows that this mobility 
spectrum analysis method may not necessarily be relied 
upon the provide information regarding peaks shapes or 
widths. Rather than converge on the true spectrum, with 
increasing iteration over the default model the algorithm 
appears to converge towards sharp narrow peaks. 

This algorithm has also been tested extensively on real 
data sets. For example, magnetoresistance measurements 
on the device "6016" presented in Ref. has been ana- 
lyzed. The structure of this device, grown using low- 
energy plasma-enhanced chemical vapor deposition. |34| 
is shown in Fig. 0] 

It features a strained 20 nm quantum well of pure Ge, 
on a relaxed virtual substrate of Sio.3Geo.7. The structure 
is doped with ten 5 spikes of boron 10 nm above the 
channel. This gives a sheet density in the channel at 4.2 K 
of 6.2 x 10 n cm- 2 with a mobility of 87 000 cm^^s" 1 . 

At 300 K there is significant parallel conduction 
from the doped region, giving a Hall mobility of 
1790 cm 2 V- 1 s- 1 at a sheet density of 1.8x 10 12 cuT 2 . 
We find a result for the channel mobility at 300 K to be 
2940 cm 2 V- 1 s- 1 with a sheet density of 5.7x 10 1:L cm _2 
which agrees closely with the quoted mobility value of 
3000 cm 2 V- 1 s- 1 . 



V. SUMMARY 

A powerful mobility spectrum analysis method has 
been presented, based on the application of Bryan's al- 
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TABLE II: Summary of results of BAMS analysis on the 2-peak synthetic data shown in Fig. [5] If both field directions are 
used then the peaks are found more accurately. This effect is strongest for the low mobility peak. 





Mi 


Til 


M2 


n 2 






[10 12 cm" 2 ] 


[cm^^s- 1 ] 


[10 12 cm" 2 ] 


Specification 


200 


10.0 


1500 


1.00 


-10T to +10T 


226 


8.8 


1482 


1.0 


OT to +10T 


377 


5.6 


1533 


0.9 



TABLE III: Summary of results of BAMS analysis on the 3-peak synthetic data shown in Fig. |H] The effects of increasing 
iteration over default models can be seen. The peak at lowest mobility is not found in the correct position because of the 
limited mobility resolution and the simple algorithm used to find the maximum. 





Mi 


ni 


^2 


n 2 


^3 


n 3 




[cn^V-V 1 ] 


[10 12 cm" 2 ] 


[cm'V-V 1 ] 


[10 12 cm" 2 ] 


[cm'V-V 1 ] 


[10 12 cm" 2 ] 


Specification 


270 


7.5 


1500 


1.00 


6000 


1.00 


1 iteration 


300 


6.5 


1421 


1.08 


5946 


1.01 


3 iterations 


300 


6.5 


1461 


1.06 


5986 


1.00 


10 iterations 


260 


7.5 


1461 


1.07 


6026 


0.99 


30 iterations 


260 


7.4 


1461 


1.08 


6026 


0.99 


100 iterations 


260 


7.5 


1502 


1.05 


6026 


0.99 




1 — I — I I I 1 1 1 1 1 — I — I I I 1 1 1 



1 iterations 
3 iterations 
1 iteration 
Actual 



A 




10- 



10" 

[i [ mW 1 ] 



10" 



FIG. 3: Mobility spectra calculated on 1000-point synthetic 
data featuring three carrier gases and 0.01% data error, with 
increasing iteration over the default model. There were 1000 
output points in the mobility spectrum (linearly spaced) be- 
tween /i = -2.0 m 2 V _1 S -1 and +2.0 m 2 V _1 s -:l but only the 
region of interest is shown. Peak positions are summarized in 
Table fTTTl On this logarithmic scale it is possible to see that 
the lowest mobility peak is only defined by a few spectrum 
points. 



gorithm to the maximum entropy method. Our method 
goes beyond a previous maximum entropy mobility spec- 
trum analvsis|l8t Il9l ] by using singular value decompo- 
sition (SVD) to reduce the dimensionality of the search 



3nm Si 

30nm Si 04 Ge 06 



10 p-type (B) 5-doping spikes 



lOnm Si 04 Ge 06 



20nm Ge QW 



40nm Sio.4Ge .6 



~10um Graded buffer Si - Si . 3 Ge 



Si(001) p-type substrate 50-100^cm 



FIG. 



n Dylow- 

energy plasma-enhanced chemical vapor deposition, [a. l34j| It 
features a strained 20 nm quantum well of Ge on a relaxed 
virtual substrate of Sio.3Geo.7- There are ten S spikes of boron 
10 nm above the channel. At 300 K there is conduction in 
both the quantum well and the doped region so mobility spec- 
trum analysis is required to extract the properties of the 2- 
dimensional hole gas in the quantum well. 



space, and consequently improve the conditioning of the 
problem. Thus, overspecification and noise in the data 
are dealt with properly. This means that a large data set 
of 100 ~ 1000 points may be used directly without a seri- 
ous penalty in terms of calculation time. User bias with 
regard to data point selection is therefore avoided. No 
other mobility spectrum analysis method can deal with 
such large data sets without a serious calculation time 
penalty. 

No previous assumptions regarding the form of the re- 
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suit are required, apart from the range over which to 
search for mobilities. Mobilities can be accurately found 
which are smaller than 1/B max . Also, since the algorithm 
works in a reduced space, there is no risk that peaks are 
lost entirely if insufficient mobility resolution is specified. 

This algorithm calculates mobility spectrum over a 
range of a values (where a controls the relative impor- 
tance of x 2 versus the entropy, that is, between a good 
fit and a physically reasonable result) and calculates the 
probability P[a] of each spectrum. Then, the algorithm 
performs an average of the spectra weighted by their cal- 
culated probabilities. The ensures that the solution is the 
best that can be reasonably deduced from the data given 
Eqs. CJ-©. In this way, the results should be relatively 
free from common artifacts such as mirror peaks, pro- 
vided that the artifacts are not introduced by systematic 
measurement errors, pi) 

However, there are a few caveats regarding mobi lity 
spectrum analysis in general. It has been suggested [lSj 
that the shapes of the peaks in a mobility spectrum of a 
non-degenerate carrier gas could give information regard- 
ing the energy-dependence of the momentum relaxation 
time and therefore the dominant scattering mechanisms. 
However, the effect on peak shape of different scattering 
mechanisms is very subtle and we believe that it is unrea- 
sonable to expect any useful results from real data. How- 
ever, if BAMS analysis is performed on data from a range 
of temperatures then the temperature dependence of the 
mobility can be examined to give information about scat- 
tering mechanisms. 

Also, whilst it should be possible to observe "harmon- 
ics" in a mobility spectrum due to anisotropy of the ef- 



fective massJ3 (for example, in the heavy- hole band of 
germanium) |22j only Ref . |24] seems to have explored this 
possibility. 

Significant intersubband scattering can also invalidate 
mobility spectrum analysis; in this case it is the intersub- 
band scattering rates rather than the transport scattering 
rates which must be considered. |35| 

One final issue to consider is that the longitudinal 
magnetoresistance in real data increases monotonically 
with field, whereas Eqs. Q-© predict that it should 
saturate. [23 This is a serious issue when analyzing real 
high magnetic field data, generic to all mobility spectrum 
techniques, and may relate to the underlying assumption 
that s(n) is independent of magnetic field. 36] Often the 
solution seems to be to assume (possibly unreasonably) 
that there are mobilities which are much smaller than 
l/-Bmaa;-l22l However, using Bryan's algorithm to balance 
the least squares fitting error with a physically reasonable 
result ensures that this mobility spectrum method should 
minimize such problems. 
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